Interplay between the single particle and collective 
features in the boson fermion model 
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We study the interplay between the single particle and fermion pair features in the boson fermion 
model, both above and below the transition temperature Tc, using the flow equation method. Upon 
lowering the temperature the single particle fermionic spectral function: (a) gradually develops a 
depletion of the low energy states (pseudogap) for T* > T > Tc and a true superconducting gap for 
T < Tc, (b) exhibits a considerable transfer of spectral weight between the incoherent background 
and the narrow coherent peak(s) signifying long-lived quasi-particle features. 

The cooperon spectral function consists of a delta function peak, centered at the renormalized 
boson energy lu — Eq, and a surrounding incoherent background which is spread over a wide energy 
range. When the temperature approaches Tc from above this peak for q = moves to a; = 0, so that 
the static pair susceptibility diverges (Thouless criterion for the broken symmetry phase transition) . 
Upon decreasing the temperature below Tc the cooperon peak becomes the collective (Goldstone) 
mode Eq oc |q| in the small momentum region and simultaneously splits off from the incoherent 
background states which are expelled to the high energy sector |lj| > 2Asc{T). We discuss the 
smooth evolution of these features upon approaching Tc from above and consider its feedback on 
the single particle spectrum where a gradual formation of damped Bogoliubov modes (above Tc) is 
observed. 



I. INTRODUCTION 

The boson fermion model (BFM) was initially invented 
for a description of the conduction band electrons cou- 
pled to the lattice vibrations in the region of intermedi- 
ate electron-phonon coupling [l| . The underlying physics 
emerges from the assumption that in the crossover region 
the electrons exist partly in form of bound pairs (hard- 
core bosons) and partly as quasi free particles (fermions). 
Many-Body correlations of such a two component sys- 
tem are induced due to charge exchange processes con- 
verting the fermion pairs into the hard-core bosons and 
vice versa. The initially very heavy bosons effectively in- 
crease their mobility and, at the critical temperature T^ 
undergo a Bose Einstein (BE) condensation while the 
fermions are simultaneously driven to a superconducting 
state. 

Upon approaching Tc from above there are several pre- 
cursor features of superfluidity/superconductivity show- 
ing up in the system. In this paper we address these 
precursor effects by means of a renormalization group 
scheme which is outlined in some detail in chapter II. We 
will show that, in the pseudogap phase (T < T* ) and 
in superconducting state (T < T^), the single- and two- 
particle properties become strongly interdependent. In 
the pseudogap phase this is seen for instance through 
a gradual destruction of the single particle states (near 
the Fermi surface) which is accompanied by a simultane- 
ous emergence of fermion pair states. Fermion pairs with 
total zero momentum show up in the normal state only 
as damped entities which propagate over a finite time 
and/or spatial scales. However, we find that there is a 
certain critical momentum qcrit{T) above which fermion 
pairs become long lived quasi-particles (which are sepa- 
rated from the incoherent background as discussed in sec- 



tion III). At the phase transition qcrit{Tc) —^ and below 
all the cooperons become good quasiparticles. The ap- 
pearance of fermion pairs above Tc affects the single par- 
ticle spectrum, leads to the formation of the Bogoliubov 
shadow branches as has been discussed by us recently 
PI . When passing through the phase transition at Tc we 
expect a smooth evolution from the pseudo- to the fully 
gaped single particle spectrum where the Bogoliubov fea- 
tures (caused by the existence of the fermion pairs) are 
present below as well as above Tc- In the superconduct- 
ing state we find that the interdependence between the 
single- and two-particle properties leads to the character- 
istic peak-dip-hump structure. Similar conclusions have 
recently been reached independently by Fieri et al. Q 
within a different theoretical approach. 

As has been widely emphasized by Uemura j^], there 
are a number of very convincing experimental indica- 
tions for precursor phenomena in the underdoped HTS 
cuprates. Whether the whole pseudogap phase can be ex- 
clusively attributed to such precursor effects is still un- 
der debate. Nevertheless, for temperatures sufficiently 
close to Tc (in the underdoped samples) the existence 
of fermion pairs, being correlated on a small spatial and 
temporal scale, were confirmed by the measurements of 
the optical conductivity in the terahertz regime J5(. On 
the other hand the static experiments measuring the 
Nernst coefficient Q gave indications for the existence of 
"moving pairs" with a certain phase slippage above Tc. 
These facts together with the peak-dip-hump structure 
found by the ARPES Q acquire here a natural explana- 
tion within precursor phenomena which are intrinsic to 
the BFM or similar scenarios, accounting for strong pair 
fluctuations. 

On a more general basis, the BFM is often believed 
to capture essential aspects of the crossover physics be- 



tween weakly coupled and strongly paired lattice elec- 
trons H, 13 ■ Various unconventional properties of the 
superconducting state have been investigated within this 
model by several groups '!(}, HAl [13, [13 ■ Several authors 
concluded according to phenomenological considerations 
[ij, [13 that the BFM can serve as an effective model for 
the description of quasi 2-dimensional strongly correlated 
cuprate superconductors jig . 

The crossover issue and the BFM turned out to be 
of particular interest also in atomic physics [171 [iM [l9l| , 
where a resonant Feshbach scattering is induced between 
the trapped alkali atoms, such as ""^K or ^Li. By ap- 
plying external magnetic fields the effective interaction 
between atoms can be varied from the weak (BCS) to 
the strong coupling (BE) limits. Under optimal condi- 
tions a resonant superconductivity is expected to arise at 
Tc ^ 0.5Tp [13, which is presently routinely observed 
in several research labs |20j. This very general scenario 
of the Feshbach resonance can be theoretically expressed 
via the BFM, as was recently shown by one of us [2ll|. 

It has been frequently stressed in the literature [2J, 
[2a . [2q| that a selfconsistent and conserving treatment of 
single-particle and pair correlations has a crucial impor- 
tance for the description of the HTS cuprates. In this pa- 
per we study the mutual interdependence between such 
single- and two-particle properties (paying special atten- 
tion to the precursor features) by extending our previous 
work ^21 based on the flow equations method [23| ■ Our 
former study focused on a diagonalization of the Hamilto- 
nian and the determination of the renormalized fermion 
and boson energies [22|. In the present paper we derive 
the Green's functions (dynamic quantities) which deter- 
mine the propagation of single fermions, single bosons 
and of fermion pairs. From these functions we obtain 
the corresponding excitation spectra. The methodologi- 
cal virtue of the flow equations method is that, besides 
treating the single- and two-particle entities on equal 
footing, it distinguishes between the contribution of long- 
lived and damped quasi-particles in the spectrum. The 
former are usually represented by delta function peaks 
with a given spectral weight while the latter are given in 
form of a broad incoherent background. 

For our study we use the following Hamiltonian [jj] 



H = ^ (£k - /i) cLcka + ^ (As - 2//) fe^6q 



k,(T 



-7^ ^(foljCq_kiCkT+h.C.) , 
k,q 



(1) 



where the operators c]^^ (ck^) refer to creation (annihi- 
lation) of fermions with the energy ek and &^ (6q) corre- 
spondingly to bosons in localized states As. The boson- 
fermion coupling v will be taken here as isotropic, al- 
though for the real HTS systems it should be used with 
a d wave prefactor [T^, [iJ, uB, ■ For simplicity we ne- 
glect here also the hard core property of bosons, which is 
justified as long as the concentration of bosons is small. 



II. THE METHOD 



Generalities 



We apply a canonical transformation S{1) in order 
to eliminate the interaction between the boson and the 
fermion subsystem. This transformation will be car- 
ried out in a continuous way (I denotes the continuous 
flow parameter) so, that the transformed Hamiltonian 
H{1) = e'^^'^-ffe"'^^'^ reduces to a manageable form for 
further analysis. The more generally known classical sin- 
gle step transformations projects out the terms which are 
linear with respect to a given perturbation. Here we de- 
mand much more stringent constraints on a transformed 
Hamiltonian going beyond such a standard perturbative 
scheme. 

The evolution of the Hamiltonian H{1) with respect to 
the varying flow parameter I is determined through the 
differential equation 



dH{l)/dl^ [t]{1),H{1)] 



(2) 



subject to the initial condition H{Q) = H. A generating 
operator is defined by ri{l) = (de^^''^ /dl) e~^^^\ 

In principle, one can transform the Hamiltonian in 
many different ways by choosing various operators ri{l) 
(or S{1)). Some particularly efficient schemes have been 
proposed by Wegner [23 and independently by Wilson 
and Glazek [23 going back to the RG approach ideas 
[23 ■ Through a continuous transformation of the Hamil- 
tonian one effectively renormalizes its coupling constants 
while keeping a given constrained structure. In other 
words, the parameters of the Hamiltonian such as the 
energies, the two-body potentials and so on are assumed 
to be Z-dependent. 

In some distinction from the RG approach one does 
not integrate out the high energy states, but instead of it 
they are renormalized in the initial part of transformation 
until / ^ (Ae)~^ [23|- Subsequently, states with small 
energy differences start to be renormalized and finally, for 
I — > 00, the transformed Hamiltonian H(oo) eventually 
reduces to a (block) diagonal structure [30l |. 



B. Implementation to BF model 

In our previous work |23 | we have derived the continu- 
ous canonical transformation for block diagonalization of 
the BFM. The transformed Hamiltonian was constrained 
to the following structure H{1) = i?o(0 + Hint{i), where 

Ho{l) - ^(£k(0-M)cLcka + 5](i?q(0-2M)6^6q 

k,(T q 

+ j;^ Yl t^k,p,q(OCkTCpiCqiCk+p-qT , (3) 

k,p,q 

H„uil) = ^=^Wk,p(0(&p+kCkiCpT+h.c.) . (4) 



k.p 



For elimina ting: Hint{l) we follow the idea proposed 
by Wegner J23 who showed that, upon using rj{l) = 
[Ho{l),Hint{l)], one obtains lim;^oo Hmt{l) = 0. In this 
case 



vii) 



-7^ Y^ ak,p(0 (&p+kCkxCp| - h.. 

k,p 



(5) 



where ak,p(0 = (ek(0 + £p(0 - -Bk+p(0) ^'k.p(0- AH l- 
dependent parameters of the Hamiltonian (|3I4|I are deter- 
mined via a set of the flow equations (16-21) given in Ref. 
[23 |. They are obtained from the operator equation 
by reducing the higher order interactions through normal 
ordering (linearization). 

Since at Z = c» the bosons are no longer hybridized 
with fermions, we essentially obtain the (semi) free sub- 
systems with renormalized effective spectra. Only the 
fermion part contains the long range Coulomb interac- 
tion C/k,p,q which in some cases can play an important 
role. For instance, in the q = p channel one obtains 
[21| a resonant type amplitude of the potential C/k.p.p(oo) 
for Ek + ffp = -Ep+k- This corresponds to the resonant 
(Feshbach) scattering between electrons when their total 
energy is equal to energy of the bound pair (hard-core bo- 
son) [13 • In the context of HTS such unusual scattering 
is also important leading to the particle-hole asymmetry 
of the low energy spectrum both, in the pseudogap and 
superconducting phases |3ll |. 

One should remember that the amplitude of the in- 
duced Coulomb potential t/k,p,q(oo) is finite for any chan- 
nel. It was estimated to be residual, of the order v"^ [2a| . 
In the following we will thus treat the transformed Hamil- 
tonian 



H{oo) = H^{oo) + H^{oo) 



(6) 



as composed of two contributions from bosons H^{oo) = 
Sq -^q^q^q '^ith their effective energy Eq = i<^q(oo) — 2^, 

and fermions approximated by H^{oo) ~ X^ko- ^kCj(.g.Ck<T 
with the effective energy 



£k = £k(oo) - fJ-+j^Yl ^^■pA°°)^p, 



(7) 



Here n^ ^ ~ (c|, a-'^p.cr): where the spin is a dummy index 
which will be kept throughout the remainder of this paper 
in order to indicate the origin of such terms. 



usual, T denotes ordering with respect to the imaginary 
time T = it. 

The computation of the thermal averages (...) = 
Tr {e^^^ ...} /Tr {e^^^} is easiest to carry out using the 
transformed Hamiltonian H{(X)) because of its (block-) 
diagonal structure. Due to the invariance of the trace 
under the unitary transformation we can write 



Tr 



{e-^^(')O(0} , (9) 



= Tr 



where 0{l) — e'^'^'-'Oe^"^'-'^. Hence, if we want to use the 
transformed Hamiltonian H(l = 00) in the Boltzmann 
factor e"''^^'-' we ought to transform the observable O 
too. For the continuous transformation this is however 
a nontrivial problem because, in order to get O(oo), one 
must analyze the whole transformation process. The evo- 
lution of the arbitrary observable 0{l) with respect to I 
must be deduced on a basis of the differential equation 



dO{l)/dl^ [?/(/), 0(0]. 



(10) 



For the Hamiltonian O = H li thus is given by (01 which 
was already discussed previously by us 22] for this model. 
In the next sections we study the /-dependence of the 
individual boson and fermion operators and of fermion 
pair operators. By looking at the limit Z — > cxj, we shall 
derive effective spectral functions 

^F,i3,pa»r(j^^^-) = -1 Im G^'^'^"'' (k, w), (11) 

IT 

where G(k,tj) = j^ dr e™G(k, r) with the following 
single particle Green's functions 

G^(q,r) = ((6q(r);6^)), (12) 

G^AKr) = ((ckT(r);c[,)), (13) 

Gf,(k,r) = ((ckT(r);c.kx)), (14) 

and the two-particle pair propagator given by 

G'^"'-(q,r) = i^^((ckx(r)cq_kT(r);<_p^c|,^)).(15) 

k,p 

We will next investigate the structure of these spectral 
functions and discuss their related physical properties. 



III. BOSONS AND COOPERONS 



C. Dynamical quantities 



A. Flow of the boson operators 



In this work we focus on determining thermal equilib- 
rium Green's functions of the form 



((Oi(T);02))=-(r,Oi(T)02) 



(8) 



where the time evolution of the operators is given by 
0{t) = e'^Oe-^", with r S (0,/?) and (3 = l/fcsT. As 



In the course of such a continuous transformation, the 
initial boson operator 6q becomes convoluted for Z > 
with the fermion pair (cooperon) operator. This can be 
seen from the / = derivative 



dbg 
dl 



[^(o),fc. 



1 



Xl"k,q-k(0)ckiCq_kT- (16) 

k 



Physically this means, that while disentangling the boson 
from fermion subsystem we obtain some new quasiparti- 
cles made out of the initial bosons and cooperons (like in 
the BCS theory where the quasi-particles are composed 
of electrons and holes). 

Guided by the structure of to H16(l it is judicious to 
choose the following superposition for the ^-dependent 
boson operator 

&q(0 = ^q(0&q +^J2 ^q,k(OckiCq_kT , (17) 

where Aq{l), Bq^k{l) are some complex functions with 
the initial condition .4q(0) = 1 and Sq,k(0) = 0. Substi- 
tuting H17|) into the flow equation (|10|l we obtain 



dl 
dl 



^E"k.q-k(0/k.q-kSq,k(0, (18) 



N 



ak,q-k(0-^q(0: 



where we introduced the shorthand notation 



^k^i 



/k,p ~ 1 - "ki ~ "-pT- 
From (|18I19|I we notice the following invariance 



l^q(Ol' 



N ^^ 

k 



|Sq,k(Ol'/k.q-k = 1 



(19) 



(20) 



(21) 



which guarantees that the commutation relations be- 
tween the /-dependent boson operators [6q(/), fe|,(/)] ~ 



^q,p 



are correctly preserved. 



The parameterization H17|) which follows from the flow 
equation (|10() for the operator bq{l) yields the following 
boson spectral function 1)11(1 



A^(q,u;) = \Aq{^)\H [uj - Eqj (22) 

+ -77 X! /k,q-k|Sq.k(oo)P^ (cj - £k ^ £q-k) • 



The first term in equation H22|) describes the coherent 
part of the boson spectral function corresponding to the 
long-lived quasi-particles with the renormalized energy 
Eq. The second contribution describes some incoherent 
background of the boson spectral function which repre- 
sents the states of a relatively short life-time. 

We solve the flow equations (|18I19|I fully selfconsis- 
tently applying the numerical procedure based on the 
Runge-Kutta algorithm. For any value of / we discretized 
the coefficients Aqil), '8q,k(0 using a mesh of 4000 
equidistant points for representing the vectors k and q in 
the Brillouin zone. Due to computational limitations we 
restrict ourselves to a bare one-dimensional tight binding 
dispersion ek(/ = 0) = — 2icos(A;a;a) and throughout this 
paper use the bandwidth D = 4t as a unit for energies 
and for the temperature. Starting from the initial value 
^q(O) = 1 the /-dependent coefficients are calculated 
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FIG. 1: The boson spectral function ^ (q = O.oj) for 
three representative temperatures T = 0.02 > T* (top), 
T — 0.007 < T* (middle) and for the ground state super- 
conducting phase T = (bottom). The BFM parameters are 
ntot = 1, V = 0.1, As = —0.6 (energies are expressed in units 
of the initial fermion bandwidth D = 1). 



via the following scheme Aq{l + SI) = .4q(/) + SI — -^i, 
where the derivative is given in equation H18() . The co- 
efficients y8q^k(0 ^^^ determined in the same way. Since 
the renormalizations of both these coefficients (as well as 
other quantities such as energies Sk{l), Eqil) and boson- 
fcrmion coupling Wk,p(0) occur at the initial steps of the 
transformation procedure, we adjust the increment 51 in 
the following way: SI = 0.01 (for ^ < 5), (5/ = 0.1 (for 
5 < ; < 10^), SI = 1.0 (for 10^ < / < 10^), and SI = 10 
(for 10'^ < / < 10^), where both / and SI are expressed in 
units D^^ . The asymptotic (fixed points) are obtained 
already around / ~ 500 but the transformation procedure 
is continued up to a good convergence i.e., I = 10^. 

Figure Q shows the results obtained numerically for 
the single particle boson spectral function A^{q,uj) in 
the long wavelength limit q = 0. We illustrate three 
distinct situations corresponding to: the normal phase 
above T* (top panel), the normal phase with the pseudo- 
gap structure present in the single particle fermion spec- 
trum (middle panel) and the superconducting state at 
T = (bottom panel). In this paper we have chosen the 
same set of parameters as previously [22, i.e. Ab = 0.4, 
ntot = I, V = 0.1 such, that the temperature at which 
the pseudogap begins to open up is roughly T* ^0.1. 

Our study of the superconducting phase (following the 
previous work |22|') is based on a three-dimensional sys- 
tem with a BCS type of approach, as far as the fermionic 
subsystems is concerned, and a BE condensation ap- 
proach for free bosons, as far as the bosonic subsystem 
is concerned. We notice, that in this phase there is a 
perfect separation of the coherent part (describing the 
long lived quasiparticles) from the incoherent part of the 
spectrum. Moreover: (a) the coherent peak is pinned at 



Lu — allowing for a macroscopic occupancy of the zero 
momentum state by a certain fraction n^g^d ^^ ^^^ ^^ 
condensed bosons, (b) the incoherent part A^^(q, w) ex- 
ists only outside a energy window (equal 2v yn^g^d ^^ 
will be explained in section IV. D). Owing to such a be- 
havior, the condensed bosons are not damped and they 
are able to establish a long range order parameter in the 
boson subsystem. On the other hand in the normal phase 
above T* the coherent and incoherent parts overlap with 
each other and consequently the boson quasiparticles are 
damped. This damping is caused by some very reduced 
remanent inter-boson interaction of the order of w*, which 
arises in this renormalization procedure 22] . 

The pseudogap phase (middle panel of Fig. ^ repre- 
sents some intermediate situation, where we notice that 
the incoherent background is partly pushed away from 
the coherent peak. Thus the zero momentum bosons 
start to emerge as better and better quasiparticles upon 
approaching Tc from above. Yet, the zero momentum 
boson state is macroscopically occupied only below Te- 
la the BFM there is a strict relation between the sin- 
gle particle boson and fermion pair excitation spectra (see 
equations H29I30|I in the next section). By inspecting Fig. 
1 (and figure 3 presented below) we conclude that zero 
momentum fermion pairs gradually emerge in the pseu- 
dogap phase (T* > T). Upon lowering the temperature, 
the surrounding incoherent background fades away and 
thus effectively leads to increase of the life-time of the 
zero momentum bosons and fermion pairs. For T < Tc, 
these entities acquire an infinite life-time. Experiments, 
sensitive to the short lived Cooper pairs, should be able to 
detect their presence above Tc. This type of a precursor 
phenomenon was indeed observed for the HTS cuprates 
using the alternating magnetic fields in the terahertz fre- 
quencies regime l^- A residual Meissner effect was seen 
there up to nearly 25 K above the transition temperature 
Tc and which is an indication that propagating fermion 
pairs exist there on a corresponding short time scale. 

In figure [5] we compare the boson spectral function 
A^{q,uj) for several momenta q = |q| at two different 
temperatures: above (column on the left) and below T* 
(column on the right). We notice that at finite momenta 
there is is a qualitative difference between these spectra. 
At temperatures T > T* the coherent peak is always cov- 
ered by the incoherent background signaling that these 
components are convoluted with each other. On the con- 
trary, at temperatures below T* , the coherent boson peak 
separates from the incoherent states. Such a splitting-off 
is very sensitive to moderate temperature changes and, 
for example at T = 0.004, occurs above a critical mo- 
mentum Qcrit — 0-01 X ir/a. This critical value decreases 
with decreasing temperature and finally qcrit(Tc) -^ 0. 
We thus see that for momenta q > qcrit{T) the coher- 
ent boson states, representing the long lived propagating 
modes are not scattered by the incoherent background. 

Our finding, that under certain conditions the boson 
peak becomes well separated from the background of in- 
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FIG. 2; Variation of the boson spectral function A (q, a;) 
with respect to changing momentum q at high temperature 
T — 0.02 > T* (panels on the l.h.s.) and low temperature T = 
0.004 (panels on the r.h.s.) corresponding to the pseudogap 
region T* > T > Tc- The lattice constant is taken as a unit 
a = l. 



coherent states ^^^^.(q, w), means that in the pseudogap 
regime {T* > T > Tc) the finite momentum bosons are 
well defined quasi-particles with infinite life-time. Since 
bosons are closely related to the fermion Cooper pairs, 
we further conclude that for q > qcrit {T) there exist infi- 
nite life-time cooperons above the transition temperature 
Tc- According to the Heisenberg principle, experiments 
measuring the pair-pair correlations restricted to spatial 
distance smaller than Ax ~ l/qcrit{T) (of the order of 
several lattice constants) should be able to detect such 
fermion pairs. And indeed, such long-lived "moving" 
fermion pairs were observed up to very high temperatures 
above Tc by applying a thermal gradient and measuring 
the Nernst coefficient in underdoped HTS materials p. 



B. Flow of the cooperon operators 

We have seen in the previous section that the boson 
operators get mixed for / > with the fermion pair op- 
erators. Now we expect that, in turn, also the latter be- 
come convoluted with 6q during such a transformation. 
By calculating the initial (1=0) derivative of the cooperon 
operators S^ = N^^ J^k CkiCq_kT we obtain 



d5q 
dl 



y^ ak,q-k(0)/k,q-k bq 



(23) 



In accordance with the previous substitution H17() we pro- 
pose the following Ansatz 



;= J2 -^k^qCOckiCq-kT + AAq(06q (24) 



5q(0 



with the initial conditions A^q.k(O) = 1 and A/'q(0) — 0. 
After substituting the expression (|24|) into (|10|l we get 



dl 
dl 



ak,q-k(OA/q(0 (25) 

^5]ak,q-k(0/k,q-kXq,k(0(26) 



N 



We recognize that the flow equations (^51 EHl for the 
unknown coefficients Mq.k(0 ^-nd Nqil) have a structure 
identical to that given by H18II9|I . There is only a differ- 
ence in the initial conditions which in this case lead to 
the following invariance 



l-^lWl' + ^E/k,q-k|A^q.k( 



In the r.h.s. of equation H27|l 
property that A^~^ I]k /k,q-k 



= \-nF. (27) 



np 



N-'Tv 



•"ka 



we made use of the 
= 1 — np, where 
denotes the total concentration of 



fermions. Equation (|27(l assures the proper statistical re- 
lation between the cooperon operators [S'q (/), 5*+ (?)] = 

7V~i X]k(l ^ '^q-kt'^q-kT ^ cj^i^ki), which can be approx- 
imated by the c- number ~ I — np. 

With the Ansatz H24() and using its hermitian conju- 
gate S't{l) we can now determine the two particle fermion 
Green's function G'P°"'(k, r) defined in IjlSI) . The corre- 
sponding spectral function l|ll|) becomes 



AP"'^(q,o.)-|AAq(oo)|25L,-^, 



(28) 



-JT X! /k.q-k|A^q,k(00)p(5 (W - £k - £q-k) 



The cooperon spectral function turns out to have a struc- 
ture similar to A^{(\,uj), expressed in H22|l . This is a 
general feature of the BFM which, in particular, im- 
plies that bosons BE condense simultaneously with the 
fermion pairs driven to superconductivity |82l |. 

We limit our quantitative discussion of the cooperon 
spectral function lP5|l by only presenting in Fig.|2|a dis- 
tribution of the occupied coherent and incoherent states 
in the long wavelength limit q = 0. Above T* almost all 
the cooperons exist as incoherent objects. Upon decreas- 
ing the temperature more and more incoherent pairs are 
converted into coherent ones. Yet, because of their over- 
lap with the incoherent background the Cooper pairs can 
propagate only on a short time scale. 

The fact that the overall structure of the cooperon 
spectral function ^^"^''(q, w) is to great extent similar to 
the single particle boson function A^(q, lo) is not surpris- 
ing. Using for instance the equation of motion technique 
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FIG. 3: The cooperons' occupancy [exp(Lj/fcflr) — 1]~ x 
A''""^(q,a;) of the zero momentum coherent (left h.s. panel) 
and incoherent (right h.s. panel) states for three represen- 
tative temperatures: T = 0.02 (above T*) and T = 0.007, 
0.004 (below T*). With lowering temperature the incoherent 
fermion pairs are converted into the coherent quasi-particles. 



for the Green's functions one can prove the following im- 
portant identity |32| 

G^(q,c.) = G^(q,c.)+ (29) 

+ t;2G^(q,c.)GP"'-(q,c.)G^(q,c.), 

where Gq{(\,uo) = [w — A^ -I- 2/.t] ^ denotes the single 
particle boson Green's function of the non-interacting 
system. Both Green's functions G^"*'"(q, w) and 
G^(q, w) have thus common poles and their spectral 
functions are correspondingly related via 



y4P"'^(q,t^) 



(c 



2m)' 



A^(q,c.) 



(30) 



For q — 0, when the single particle boson Green's func- 
tion develops a pole at w = 0, the two particle cooperon 
propagator is characterized by the same pole - although 
with a different weight. This automatically causes a 
divergence of the static pair susceptibility and via the 
Thouless criterion leads to the phase transition into a 
superconducting state. 

We want to stress that all our conclusions concerning 
the evolution of the boson spectral function A^ (q, lS) are 
also valid for ^^'"''(q, w). In particular, we emphasize, 
that in the pseudogap region T* > T > Tc there ex- 
ist fermion Cooper pairs. Fermion pairs with total zero 
momentum are convoluted with incoherent background 
states and therefore their life time is finite. They can 
be detected experimentally only via some short time im- 
pulses. On the other hand, long-lived fermion pairs can 
safely exist if their total momentum exceeds the critical 
value q > QcritiT). In this case correlations between the 
fermion pairs can possibly extend over the spatial dis- 
tance up to Ax ~ ^/qcrit- A possibility to observe corre- 
lations between the fermion pairs above Tc over a short 



spatial and temporal scale had been previously suggested 
by Tchernyshyov p^ . 



IV. FERMIONS 

A. Flow of the single particle fermion operators 

Some aspects of the effective single particle fermion ex- 
citations were already discussed in our recent Letter ^. 
We focused our discussion there on the single particle 
fermion spectrum for a narrow region around the Fermi 
surface k ~ kp. We have shown that the pseudogap 
is accompanied by an appearance of the Bogoliubov-like 
branches. The shadow part of these Bogoliubov modes is 
visible above T^ as a broad structure which narrows upon 
approaching Tc- Below Tc both branches become in- 
finitely narrow, signaling that bogoliubons become then 
the true, long-lived quasi-particles. This result is here- 
with confirmed from our present direct study of the boson 
and fermion pair spectra, illustrated in figures 1-3. 

In this section we would like to present a detailed 
derivation of the diagonal and off-diagonal parts (in 
a Nambu spinor representation) of the single particle 
fermion Green's function. We will discuss the diagonal 
and off-diagonal parts studying their structure over the 
whole Brillouin zone above and below r^. 

Let us briefly recapitulate the main properties of the 
fermion operators cj^.^ and Cko- resulting from the contin- 
uous canonical transformation. At the initial step (i.e. 
for / = 0) their derivatives read 

rfCkcr ,„^ bo + 



-^ E "q-k.k(0)foq<_k.-. (31) 



q#0 

where the negative sign corresponds to cr =t and the 
positive sign to a —[. The first term on the r.h.s. of 
(|^ shows that the fermion particles get mixed with the 
fermion holes. The last term in l|31(l corresponds to scat- 
tering of fermions on bosons with finite momentum. 

As a consequence of (|31ll we imposed the following 
parametrization of the /-dependent fermion operators 

CkT(0 = ^k(0 CkT + 7^k(0 c^_ki (32) 

+ 7^7^ Pk,q(0 ^^]jCq+kT + rk,q(0 foqc|j_ki 

clki(0 = - T^lil) CkT + K{1) cLki (33) 

~ 7^ H [''k.q(0 feqCq+kT " Pk,q(0 ^qCq„ki " 
^ q#0 

In this case, the initial conditions read 7^k(0) = 1, 
7^k(0) = and Pk,q(0) = 0, rk,q( 0) -^ 0. As we 
have stated previously ^, equations (|32I33|I generalize 
the standard Bogoliubov- Valatin transformation 33] in 



a twofold way: (a) the initial particle and hole operators 
are transformed into cooperons (bogoliubons) through a 
continuous transformation, (b) the scattering of finite 
momentum Cooper pairs is additionally taken into ac- 
count via terms containing the coefficients Pk,q(0 ^^d 

''k,q(0- 

/^From ()10|l one can derive the set of differential equa- 
tions for a determination of the /-dependent coefficients 
appearing in H32I33(I . Because of their importance, for 
our further discussion we repeat them here again [2| 



dl 



dUkil) 
dl 



rfpk,q(0 
dl 

drk,ci{l) 



.id a-k.k(0 ^k(0: 



(34) 



J^ Yl "q-k,k(0 ('^q + "q-ki) ''k,q(0 



q#0 



ak,-k(/) ^k(0, 



(35) 



-T7 Y^ a-k,q+k(/) (nq + n^+kt) Pk,q(/) 



dl 



N 

q#0 

a-k,q+k(0 T^k{l), 
- ak,q-k(0^k(0- 



(36) 

(37) 



7^ 

cond 



bosons and n^ = (&l&q) is the distribution of finite 
momentum q 7^ bosons. Equations H34I37|I prop- 
erly preserve the fermionic anti-commutation relations 

{ckT(0>CpT(0^} = <^k,p and {c\^a{l),Cpa'{l)} = 0. 

The single particle fermion Green's function can be 
easily obtained in the limit / — > 00 when the fermions are 
no longer coupled to the boson subsystem Wk,p(oo) = 0. 
We have shown previously Q] that the diagonal part reads 

A^AKu^) (38) 

= \V^{^)\''5 {uj - ek) + |7^k(cx3)|2,5 {uj + e_k) 

7 XI ('^q + "q+kt) bk,q(00)|^(5(w + Ec^- Sq+k) 



N 



q#0 



TV /-^ V q 



^q-ki) kk.q(00)|2(5(w -E^ + eq_k), 



q#0 



and that it correctly satisfies the sum rule 
^ dLoA^{k,uj) — 1. The spectral function l|38|) con- 
sists of the delta function peaks which represent the 
long-lived quasi-particles and a remaining incoherent 
background A^j„^(k, cj) of the damped fermion states. 
It should be mentioned that a similar result was obtained 
by M. Fisher and coworkers 34] for the 2 dimensional 
Hubbard model with a fractionalized structure. Those 
authors considered the spinon and chargon degrees of 
freedom coupled to a Z2 gauge field. At low temper- 
atures the spinons and chargons are confined. Their 
deconfinement becomes possible at finite temperatures 
by overcoming the gap of "vison" excitations. From 
an analysis of the low energy excitations those authors 
derived the effective spectral function for physical 



electrons which took exactly the same form as our result 
(|38|l . Apart from a common structure of the spectral 
functions, the remaining discussion and interpretation 
for both models is different. 

The spectral function for the off-diagonal part of the 
fermion Green's function ((ckt('''); c_kx)) becomes 

Afd(k, w) = 7'k(oo)7^k(oo) [S (uj + Ek) -5{lj^ Sk)] 



'iV ^ Pk,q(00)rk^q(00) |(5(cJ + Eq - Eq+k) 



(39) 



q#0 



X ['^q + "q+kt] - ^(^ -E^ + £q-k) [t 



'q-kiJ 



which conserves the sum rule /da;A^^(k, cj) — 0. This 
function consists also of the delta peaks and some frac- 
tion distributed over a wide energy region. In what fol- 
lows below we shall study the properties of the diagonal 
and off-diagonal spectral functions in various tempera- 
ture regions. 



B. The normal phase spectrum 




FIG. 4: Dispersion of the effective single particle fermion en- 
ergy ik for several temperatures (upper panel) and spectral 
weights of the coherent and incoherent parts for T = 0.02 and 
0.004 (bottom panel). Momentum is expressed in units 1/a. 



Above the critical temperature Tc there is no boson 
condensate nf^^^^ = and in such a case the flow equa- 
tions (|34I37I) become decoupled from (|35I36|I . Conse- 
quently the coefficients TZk{l) and Pk.q(0 '^^ ^ot change 
their initial zero values, i.e. 



7^k(0 = 0, 



Pk,q(0 = 0. 



(40) 



This property (|40() has a strong influence on both fermion 
spectral functions. The off-diagonal part identically van- 
ishes A^^CkjUj) — what implies that above Tc there 
exists no long range order parameter in the fermion sys- 
tem. 

The spectral function of the diagonal part is described 
by 



A^(k,c.)-|7'k(oo)|25(a;-£k) 
1 v^ 



(41) 



— / ("n 

TV /-^ V q 

q#0 



"q-ki) l?'k,q(00)p,5(cj -Eq + £q_k) 



We thus have just one branch of long-lived fermion states 
described by the dispersion Ck and which can be obtained 
from the renormalization scheme discussed by us in Ref. 
[2^ . For a clear understanding of the resulting low energy 
physics we show the temperature evolution of this quan- 
tity in figure ^ together with the corresponding spectral 
weight |7'k(oo)p. Below a certain characteristic temper- 
ature T* (in this case T* ^ 0.1) we observe a tendency to 
form the pseudogap. Simultaneously there occurs a par- 
tial transfer of the spectral weight from the quasi-particle 
peak to the incoherent background states (see the bottom 
panel of Fig. 0}. 

Such a transfer of the spectral weight is responsible for 
the appearance of the Bogoliubov shadow band, as shown 
from the selfconsistent numerical calculations in Ref. 13 . 



We will prove this result here analytically. When the 
temperature approaches Tc from above, bosons tend to 
occupy macroscopically the lowest lying states of energies 
Eq ~ 0. These states are spread over the momentum 
region |q| < A with the cutoff A being a fraction of the 
inverse lattice spacing a~^. Since the boson occupancy of 
low momenta is much larger than the fermion population 

F 
"q-k 



we can simplify H41|l to the following form 



A^iKco) c |7'k(^)|2 5(a;-£k) 



B 



|r-k.q^o(cx3)P Y. -^S{u; + e.^-Eq) 
|q|<A 



+ ^d,riqtdV^^ ^) ^ 



(42) 



where A 



d,rigid\'^^'^l — ^' ZL^|q|>A"-q-kikk,q(oo)| 

q_k — -Eq) represents a rigid background which is 



weakly sensitive to varying the temperature. Integrating 
over the small momentum boson states (whose energies 
are negligible E^^ ~ 0) we obtain a partly broadened peak 
with its maximum occurring at w = — £k- This additional 
branch of fermion excitations can well be fitted by the 
Lorentzian shape 



E 

|q|<A 



-f"^" 



£-k 



-Eq) 



n 

TT 



r,. 



(c. 



£-k) 



n 



(43) 



with n^ being the concentration of bosons. For T < Tc 
(|43|l shrinks to F ^ 0. Similar Bogoliubov shadow bands 
were also indicated for the normal phase of the f^ < 
Hubbard model using the Quantum Monte Carlo studies 
35] and the conserving diagrammatic treatment of Trem- 
blay and coworkers fSQI- In our case the Bogoliubov- like 



A;(k,co) 



QJ-K 



k=0.6ii 



k=0.57i 




FIG. 5; Diagonal part of the single particle fermion spectral 
function 1411 1 for a large part of the Brillouin zone at T = 
0.007. Vertical bars represent the coherent part A^^^^y^{\i.,uj) 
while the shaded area corresponds to A^ i^^iV., uj) . Besides 
only a close vicinity of the Fermi surface (consult 14311 or 
see the figure 1 in Q) this spectral function is rather weakly 
affected by temperature. 



branches are characterized by the pseudogap dispersion 
±ek shown here in figure 4. Above Tc the shadow branch 
(corresponding to — sign) is broad as discussed in our 
Letter 12 . 

Far away from the Fermi surface the renormaUzed en- 
ergies are nearly equal to the bare values Sk = e\^ — jjl. 
The long-lived quasi-particle contains however only part 
of the spectral weight, i.e., |7-'k(oo)P < 1. In figure [S] 
we illustrate the coherent and incoherent contributions 
of the single particle fermion spectral function 1)41(1 in 
a large part of the Brillouin zone. The damped fermion 
states are spread over the large energy regime lo, they ex- 
ist even for momenta quite distant from kp. This result 
agrees well with the previous studies based on selfconsis- 
tent perturbative theory (see the third reference of [l^). 
The presence of the substantial incoherent background 
states might possibly be related to the experimental sig- 
nal observed in the ARPES measurements for the HTS 
compounds '7'|. 



C. The long-lived quasiparticles below Tc 



B 



perconducting phase the effective dispersion is thus char- 
acterized by the true gap l^sciT) 

For a finite condensate fraction n 
equations if. mi. 



V\J'r^cond{T)- 

^ jj, we can see from 
that the coefficients 7?.k and Pk.q be- 
come finite too. The long-lived fermion excitations are 
then represented by the two sharp Bogoliubov branches 



^d,coh(k,w) 



|7^k(cx3)|2,5(c^- 






(44) 



Similarly, the coherent part of the off-diagonal spectral 
function H39|) reads 



^fd,co/i(k,w) 



= 7'k(oo)7?.k(oo) 

X [5 (w + £k) - (5 (tj - £k)] . (45) 



Equations H44I45|) resemble the BCS-type results which 
usually app ear in the mean field studies of this model 
[lOl llll lia | . Let us stress that in distinction to these 
we obtain here the total spectral weight engaged in the 
coherent parts to be [21 



iT'kMp + |7^k(cx3)|2 < 1. 



(46) 



The rest of the weight is redistributed over the damped 
fermionic states. This will be discussed separately in the 
next section. 

In order to show the correspondence of our present 
analysis to the previous results for this BF model we 
prove in Appendix A, that upon neglecting the incoher- 
ent background states of 
usual BCS coherence factors 



one exactly obtains the 



l^k((X3)|^ 



|gk - Ml 

CMF 



1 - |7^k((^)|^(47) 



T'kloo) 7?.k(oo) = - 



cond 



2^r 



(48) 



where ^^^^ = \/{eu — m)^ + ^L- By neglecting the in- 
coherent parts, the flow equation procedure reproduces 
exactly the standard mean field equations [lOlllll ll ^ 



(clcT^kT) = 2 
(c-kiCkt) = - 



1 7TT7T- tgh 



cMF 
?k 



fMF 



2fcsT 



tgh 



tMF 



2^MF ---y2fcBT 



(49) 
(50) 



We shall generalize these equations (|49l I50|) in the next 
section, taking into account the contribution from the 
damped fermion states. 



Below the critical temperature T^ there exists a finite 
amount of the BE condensed bosons n^^^ ^ 0. This 
has important effects on the fermionic spectrum. U sing 
the flow equation method we have shown analytically |22l | 

that Ek 



9 B 



sgn(ek ~ fJ.)Jiek~ fJ.) 



+ ^■'Kond- In the su- 



D. Damped quasi-particles below^ Tc 

The finite life-time (damped) states do participate in 
the fermionic excitation spectrum both above T^, as well 
as below of it. Nevertheless their presence does not spoil 
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a long range coherent behavior between the fermion pairs 
which is necessary for superconductivity to occur. In 
chapter III we saw that boson as well as the fermion pair 
spectra are characterized below Tc by long-lived collective 
modes with i^q ex |q| 22^ which becomes separated from 
the surrounding background of incoherent states by the 
gap 2Asc- We show below that in the single particle 
excitation spectrum the incoherent states can exist only 
outside the gap A^c around the Fermi energy. 

The contribution from the damped fermionic states is 
described by the following part of the spectral function 

(EH 



q#0 



+ « + <-ki) |rk,q(00)|25(w - ^q + eq_k) 



(51) 



and similarly in case of the off-diagonal function (|5^ 

^fd.»«c(k,w) = - ^ Xl^k,q(00)rk,q(00) 



q#0 



("q + '^q+kt) ^{"^ + ^q - ^q+k) 
(n^ + »^q-ki) K^ ~ E^ + £q-k) 



(52) 



It is instructive to consider first the incoherent spectral 
functions H51I52(I for the ground state T = 0. Finite 
momenta bosons have Eq^o > 0, such that all bosons 
are condensed at T = and we can put n^ — for any 
q 7^ 0. On the other hand, fermions occupy only the 
states below the Fermi surface, i.e. n^^ = 9{—S]i)n^^. 
Therefore, the first term on the r.h.s. of (|51|l gives rise 
to the appearance of the incoherent fermion states at 
negative energies. It can be written as 



^d,mc(k,w<0) 

JJ Yl ^*q+kTbk,q(0o)p5 
eq+k<0 



(53) 



OJ - 



(^q + |eq+k| 



Similarly, the last term in (|51|l . which corresponds to 
incoherent states at positive energies, becomes 



<,nc(k,^>0) 

^ Yl "q-^kikk,q(00)p5 

eq-k<0 



(54) 



LO- E, 



|£q-k 



The off-diagonal spectral function (|52|l can be derived in 
the same manner. From inspecting their structure we 
notice that: 

• no damped fermion states are allowed to occur 
within the superconducting energy gap window 
\uj\ < Asc because |eq±k| > A^c, 




FIG. 6: Schematic view of the single particle fermion spec- 
tral function yl^j {]i.F,uj < 0) in the ground state of the boson 
fermion model. Energy uj is expressed in units of the super- 
conducting gap which is roughly scaled via Asc — 'v^/nfond- 



• even outside the superconducting gap, in a close 
vicinity of the coherence peaks H44() there is a par- 
tial suppression of the damped fermion states due 
to finite Eq in equations H53I54|) . As shown in our 
earlier studies [l^, |2^ the width of such a boson 
band is rather small and of the order v^ /D. 

The incoherent fermion states can thus be formed only 
outside the superconducting gap |w| > Asc and yet they 
are strongly suppressed over some small energy region (of 
the order v'^/D from the coherence peaks). 

To summarize our results of the last two sections , we 
conclude that the single particle fermion spectral func- 
tion H38I) is characterized in the superconducting phase 
by: (a) the presence of narrow quasi-particle peaks at 
Lu = ±\/(ek — A*)^ + ^sc: (t>) the partial suppression of 
the fermionic states very close to the coherence peaks, 
giving rise to the appearance of the dip-like structure, 
(c) the existence of a broad incoherent background spec- 
trum with its flat maximum (hump) situated fairly away 
from the Fermi energy. Upon increasing the temperature 
we expect the dip-like structure to be gradually filled in 
via absorbing part of the spectral weight from the coher- 
ent Bogoliubov peaks. Our expectation is motivated here 
by the fact that the spectral weights |7'k(oo)p, |7?.k(oo)p 
of the coherent peaks are very sensitive to temperature 



because of n-^ondi'E) appearing in their flow equations 

Our theoretical prediction concerning the single parti- 
cle spectrum of the superconducting phase is in qual- 
itative agreement with the experimental data known 
for the high Tc underdoped cuprates. The direct pho- 
toemission spectroscopy probes the occupied spectrum 
A^CkjLu) X [1 + exp{uj/kBT)]~ . Photoemission mea- 
surements indeed revealed the appearance of such peak 
- dip - hump structure 7] . This issue has been so far in- 
terpreted in terms of the boson resonant mode to which 
the single particle excitations are coupled. Several candi- 
dates have been proposed in the literature for such boson 
modes [33, |3a lly ■ The BF model ^ is an alternative 
scenario, where bosons represent nearly localized electron 
pairs. 
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Concluding this discussion on the coherent and inco- 
herent single particle spectra, let us stress that the expec- 
tation values for the particle distribution n^^ at arbitrary 
temperature can be obtained by integrating over the ex- 
pression 138(1 after having been multiplied by the Fermi 

Dirac function f{x) = (l + e^/'^'^-'") 



I.e. 






dL0A^{k,L0)f{uj). 



(55) 



Similarly, the order parameter determined via l|39[l is 
given by 



-kiCkt) = 7'k(^)7^k(oo) [1 - 2/(£k) 



+ JJ Xl^k,q(00)rk,q(00) [( 
q#0 

X /(eq+k 



"■q-t-kt) 



(56) 



E„ 



q-ki. 



/(-eq_k + ^q 



which generalizes the mean field equation (|50|l . Unfortu- 
nately, at the present stage we are not able to solve nu- 
merically these equations l|^l5S|l for some realistic d > 2 
system when Tc > 0. We will try to address this problem 
in the future by applying some approximate treatment. 



V. CONCLUSIONS 

In the present paper we studied the mutual relations 
between the single particle and the pair excitation spec- 
tra within the BF model QJ . We investigated their inter- 
dependence occurring in the superconducting phase be- 
low the critical temperature, and in the pseudogap region 
above Tc. Many-Body effects were treated by means of 
the continuous canonical transformation 23'] originating 
from a general framework of the renormalization group 
technique 28] . 

We have earlier reported [23 that, upon approaching 
Tc from above, there is a partial suppression of the single 
particle fermion states near the Fermi energy (pseudo- 
gap). Here, we supplement this picture by showing that 
the pseudogap feature is accompanied by a subsequent 
emergence of the fermion pair properties. Such pairs can 
show up in the pseudogap region as long lived entities 
provided that their total momentum is finite and larger 
than a certain qcrit (T) . For T ^^ Tc the critical momen- 
tum is Qcrit = 0, therefore all the fermion pairs become 
good quasi-particles. The zero momentum fermion pairs 
exist above Tc only as the damped objects because of 
their overlap with the incoherent background. However, 
upon lowering the temperature such an overlap gradu- 
ally diminishes which effectively leads to increase their 
life-time. 

Although above Tc there exist no off-diagonal long 
range order (ODLRO) we predict nevertheless the possi- 
bility to observe some ordering on a restricted spatial and 
temporal scale. Experiments using terahertz frequencies 
[5j did indeed confirm the existence of preformed pairs in 



the underdoped HTS materials up to 25 degrees above 
Tc- We moreover suspect, that "moving" fermion pairs 
(on the basis of our study expected to have the infinite 
life-time) have been observed in the measurements of the 
Nernst coefficient p|. 

As far as the single particle excitations are concerned 
we provided here the analytical arguments for the ap- 
pearance of Bogoliubov-type bands in the pseudogap 
phase. Upon approaching Tc from above the shadow 
branch absorbs more an more spectral weight and si- 
multaneously narrows as was previously indicated by us 
from the selfconsistent numerical study 2]. Below Tcj the 
shadow branch shrinks to the usual delta function peak 
and marks the appearance of infinite life-time cooper- 
ons. In distinction to the conventional BCS supercon- 
ductors we find that the quasi-particle peaks occurring 
at w = ±-\/(ek — m)^ + ^L become sfightly separated 
from the rest of the spectrum which is present in a form 
of the incoherent background. This effectively leads to 
the formation of the intriguing peak - dip - hump struc- 
ture which has been well documented experimentally by 
ARPES measurements. In our scenario such a structure 
is a consequence of the pair correlations (see also the 
similar conclusion in the Ref. Q)- 

In a future study we plan take into consideration an 
anisotropic 2-dimensional version of this model which 
is more realistic for describing the HTS cuprates. The 
boson-fermion coupling v should then be appropriately 
factorized by the d-wave form factor. Another important 
aspect which has not been fully explored so far concerns 
the doping dependence of the pseudogap and other re- 
lated precursor features discussed above. Such a study 
requires a proper treatment of the hard-core nature of 
the local electron pairs, which is an issue which is rather 
nontrivial. Preliminary results were so far obtained using 
a perturbative approach |40J and from exact diagonaliza- 
tion studies 1411. 



VI. ACKNOWLEDGMENT 

We kindly acknowledge instructive discussions with 
Prof. F. Wegner on the flow equation procedure. T.D. 
moreover acknowledges a partial support from the Pol- 
ish Committee of Scientific Research under the grant No. 
2P03B06225. 



APPENDIX 



Neglecting the incoherent background states of the sin- 
gle particle functions ^^(k,a;) and ^^^(k, w) is equiva- 
lent to the assumption that Pk,q(0 — ^^^ ^k.q(0 ~ 0. 
In such case the fiow equations (|34l I35|l simplify to 

d7'k(0 



dl 

dUkil) 

dl 



,rf a-k,k(0 7^k(0 



(A.l) 



cond 



ak, 



,(0 7'k(0- (A.2) 



12 



We will now assume that both functions Vu{l) and TZ\^{1) 
are real (this requirement does not restrict the generality 
of our considerations). We rewrite IjA.ip as 



rf^k(0 
i?k(0 



r^cond a-kM{l)dl 



(A.3) 



f-/ — oo 



and integrate both sides of ljA.3|) in the limits Jj^„ 
Using the sum rule (10) of Ref. |3] we can substitute 
7^k(0 = ^/l — ^k(0- Upon integration we get for the 
l.h.s. 



/ — OO 



dV^il) 



1=0 v/1 - (^k(O)' 



= - arcos[Pk(oo)] (A.4) 



because arcos ['Pk(O)] — 0. 

Integration of the r.h.s. of (|A.3p requires the knowledge 
of Q;_k.k(0 ^^^ the superconducting phase. First, let us 
notice that from the general definition of this parame- 
ter we have a-k,k(0 = 2\eu(l) — fi] w-k,k(0- Using the 
equation (43) of the Ref. [221 we can further write 



^-k,k 



{l)dl 



c?«-k,k(0 
4 [ek(0 - /i] ■ 



(A.5) 



In the superconducting phase we can additionally make 
use of the invariance shown in the equation (51) of the 
Ref. m 



£k(o - M - ±^ji(rr - -fo„.«^k,k(o (A.6) 

where ± — sgn(ek — m) ■ By substituting (|A.5|I and IjA.Gp 
into the r.h.s. of equation ljA.41) we obtain 



dl a- 



k,k 



(0 



,B 



/— oo 



cond I " d f_k k(v 

/=0 



ek(0 - M 



Tl 



±- 



d V- 



k,k 



/(^MF)2 



'^cond^-k.k 



arcos 



cMF 
Sk 



(A.7) 



In the second line of (jA.Tp we applied the initial condition 
i'-k,k(0) = V and also the final result f-k,k(oo) = when 
changing the integration variable from dl to dw_k,k(0- 
By comparing the results (|A.4|I and l|A.7|) we obtain 



=F 2 arcos [7'k(oo)] = - 



'"^cond 



fMF 
?k 



(A.8) 



and by taking the cosine function on both sides we get 
27'k(oo) - 1 = |ek - Ml/^k^^ which leads to the usual 
BCS coherence factors 



^kl 



oo) = i 



1 



|ek"/^| 



fMF 
Sk 



1 



7^2(oo). (A.9) 



In this way we also know that the magnitude of the prod- 
uct is |'Pk(oo)7?.k(cxD)| = 



adl ^^\i ■ 



Hence, assum- 



ing that 7-'k(0 is positive (in particular also for / — cxj) 
then on a basis of the flow equation ljA.2(l we conclude 
that 7?.k(0 = — sgn(£k — A^)|'7^k(OI ^'^^ thus finally 



Pk(oo) 7?.k(oo) 



cond 



2^MF 



(A.IO) 



This product (jA.lOp enters the off-diagonal Green's func- 
tion and in consequence yields the equation (|5()|l for the 
order parameter. 
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